Minimum variational principles for time-harmonic waves 
in a dissipative medium and associated variational 
principles of Hashin-Shtrikman type 

Graeme W. Milton 

Department of Mathematics, University of Utah, Salt Lake City UT 84112, USA 

John R. Willis 

University of Cambridge, 
Department of Applied Mathematics and Theoretical Physics, 
Wilberforce Road, Cambridge CB3 OWA, U.K. 

Abstract 

Minimization variational principles for linear elasto dynamic, acoustic, or electro- 
magnetic time-harmonic waves in dissipative media were obtained by Milton, Seppecher 
and Bouchitte generalizing the quasistatic variational principles of Cherkaev and Gib- 
iansky. Here a further generalization is made to allow for a much wider variety of 
boundary conditions, and in particular Dirichlet and Neumann boundary conditions. 
In addition minimization or maximization principles of the Hashin-Shtrikman type, 
incorporating "polarization fields", are developed. 
Keywords: waves, variational principles, dissipative media 



1 Introduction 



Stationary variational principles for wave equations have been the subject of much attention 



[see, for example, Oden and Reddy (1983), Lazzari and Nibbi (2000) and Altay and Dokmeci (2004), 
and references therein]. Among these the variational principles of Gurtin (1964), Tao (1966) 
and Willis (119811 119841) are minimizing variational principles in the Laplace domain but not 
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in the frequency domain. In the frequency domain in a dissipative medium they correspond 
to saddle point variational principles (IBorcea 1999} Milton, Seppecher, and Bouchitte 20 09). 

Saddle point variational principles for time-harmonic waves in dissipative media in the 
quasistatic limit, where the wavelengths and attenuation lengths are large compared to 
the size of the body, were derived by Cherkaev and Gibiansky (1994). They realized that 
by making a partial Legendre transform one could convert them to minimization princi- 
ples. Equivalently, one could rewrite the complex equations in a form involving the real 
and imaginary parts of the fields with a real positive definite matrix entering the consi- 
titutive law. From this reformulation, and from the differential constraints on the fields, 
one immediately obtains the minimization principles. Their variational principles were em- 
ployed to give a simple derivation (IMilton 19 90) of existing bounds on the effective com- 
plex electrical permittivity of lossy multiphase composites. [These bounds, which gener- 
alized two-phase bounds of Milton (1981) and Bergman (1982), were first conjectured by 



Golden and Papanicolaou (1985) and Golden (1986) and subsequently proved by Bergman (1986), 



Milton (1987), and Milton and Golden (1990).] Their variational principle also provided en- 



tirely new bounds on the complex bulk and shear moduli of viscoelastic two-phase composites 



(Gibiansky and Milton 1993 Gibiansky and Lakes 1993 Milton and Berryman 1997; Gibiansky and Lakes 1! 



Gibiansky, Milton, and Berryman 1999), which in turn led to bounds on the complex ther- 



mal expansion coefficient of viscoelastic two-phase composites (|Berryman 2009). 



Milton, Seppecher, and Bouchitte (2009) realized that the technique of Cherkaev and Gib- 



iansky could be extended to obtain minimization variational principles for time harmonic 
waves in a body composed of dissipative material. In particular they obtained minimization 
variational principles for acoustics, elastodynamics and electromagnetism. [In passing we 
mention that the final paragraph in that paper should be disregarded as the inequality (6.8) 
does not provide any constraints on the moduli inside the body as the difference between 
the right and left hand sides can be expressed as the integral of a square.] An unappealing 
feature of these variational principles is that they require boundary conditions on the fields 
which are difficult to physically impose, such as fixing the real part of the displacement to- 
gether with the real part of the traction (for elastodynamics) or the real part of the tangential 
component of E together with the imaginary part of H (for electromagnetism). 

Here, by including what amounts to additional surface terms, we derive more general varia- 
tional principles which allow for all sorts of boundary conditions. These include the standard 
boundary conditions where both real and imaginary parts of the displacement, traction, tan- 
gential component of E, or tangential component of H are specified. Like the variational 
principles of Milton, Seppecher and Bouchitte, these variational principles should be useful 
for tomography where one seeks information about the moduli inside a body from a series 
of measurements of the fields at the surface of the body from various different boundary 
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conditions. 



We remark that the formulation in terms of minimization principles enables one to use the 
conjugate gradient algorithm for the numerical solution using finite elements. This has been 



explored by Richins and Dobson (2010) for the Helmholtz equation, who also obtain error 
bounds on the difference between the numerical solution and the exact solution as a function 
of the grid size. 

Another completely different type of variational principle was derived by Hashin and Shtrik- 
man (1962b], I1962a[ I1962cj I1963() to bound the effective moduli of composites. This type of 



principle was extended to wave equations as stationary variational principles by Ben-Amoz (1966) 
and Willis (I1981[ 11984}) . Here we derive Hashin-Shtrikman type minimization principles for 
time harmonic waves in dissipative elastodynamic media. The extension of these principles 
to acoustics and electrodynamics exactly parallels the treatment given for elastodynamics 
and is therefore omitted. 



2 Preliminaries for elastodynamics 

The propagation of stress waves in any body is governed by the equation of motion 

diva + f = p (2.1) 

where a denotes the stress tensor, p is the momentum density (a vector) and / is body force 
per unit volume. The body is characterized by its constitutive relations which will be taken 
to have the linear form 

a = C * e, p — p*v (2.2) 
where the strain tensor e and velocity vector v are related to displacement u by 

e = i{Vm+ (Vm) t }, v = u (2.3) 

The symbol * denotes convolution with respect to time and thus the body is taken to be 
viscoelastic. Although it is usually the case that the mass density p is a real scalar, the 
formulation is not complicated if p is taken to be a tensor-valued convolution operator and 
therefore this is assumed. 

The case of time-harmonic disturbances will be assumed, in which case all fields become 
functions of the spatial variable x only, multiplied by the factor e lult which will be suppressed 
in the sequel. In this case, the superposed dot denoting time derivative is replaced by icu 
and the convolution * is replaced by ordinary multiplication. If the body occupies a domain 
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Q and tractions T = a ■ n are applied at its boundary dQ, then the rate of working of the 
surface tractions and body force is 



W = I Re(Toe^) • Re(icoue iwt ) dS + [ Re(fe iuJt ) ■ Re{iujue iwt ) dx 
Jan Jn 

= [ {Re(ae^) : Re(iuee iuJt ) + Re(iupe iuJt ) ■ Re(icoue iult )} dx, (2.4) 
Jn 

having employed the divergence theorem and the equation of motion (12. ip . It follows, em- 
ploying the time-reduced version of the constitutive relations (12. 2p . that the mean rate of 
working of the applied forces is 

W = i w / {e'C'e' + e"C"e" - u 2 (u'p"u' + u"p"u")}dx, (2.5) 



n 



where C = C' + iC", e = e' + ie" etc. Also, the relevant contractions (: or •) are left implicit. 
It follows, assuming that u> > 0, that positive-defmiteness of the quadratic forms e'C'e' and 
u'(—p")u' guarantees positive energy dissipation. This will be assumed henceforth, though 
relaxation of these conditions (such as p" = 0) will be discussed after the main development 
is complete. 

Now consider the time-harmonic problem, as described above, but with more general bound- 
ary conditions. Three compatible boundary conditions are given at each point of dQ. These 
could be for instance: all three components of displacement; all three components of trac- 
tion; two components of traction and the orthogonal component of displacement (e.g. normal 
component of displacement and shear components of traction). 

The problem with any such conditions can be described by the stationary principle of 



Willis (1984), whose time-reduced form is 

5 / {\eCe — a e — \uj 2 upu — iup u} dx = 0, (2.6) 
Jn 

the variation being taken over displacement fields u that satisfy whatever conditions are 
prescribed for displacements on dQ. The fields <7o and po are not unique but are any fields 
that satisfy the equation of motion (12.11) . together with any conditions that are prescribed 
for the tractions on dQ. To confirm this claim, note that the variation is 

{6e(a — a ) —iu5u(p — po)} dx = / 5u(a — o~o)n dS — 5u[div(a — a ) — icu(p — po)} dx. 

Jan Jn 

(2.7) 
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Requiring the volume integral to vanish for all 5u implies the time-reduced version of the 
equation of motion (12. ip . The surface integral picks out only those components of surface 
traction for which the corresponding component of displacement is not satisfied, and hence 
the boundary conditions on traction are enforced. 

There is a dual stationary principle, that 

8 / {i<7C~V — e <7 + \pp~ x p ~ iujpuo} dx — 0, (2.8) 
Jn 

where Uq is any displacement field that satisfies any given boundary conditions on displace- 
ment. The variation is taken over fields a, p that satisfy the equation of motion and any 
given traction boundary conditions. 



3 Minimum variational principles for elastodynamics 

First, the stationary principle ( 12 .6[) is spelled out explicitly: 

5 [ {i(e'CV - e"C'e" - 2e"C"e') - e'a' + e"< 
Jn 

- lu\u'p'u' - u"p'u" - 2u"p'p') + u( P y + p' u")} 
+i[\{e'C"e' - e"C"e" + 2e"C'e') - e"a' - eV£ 

- \u 2 {u'p"u' - u"p"u" + 2u"p'v!) - lo( P ' u' - p'ou")]}dx = 0. (3.1) 

Note that (12. 6p remains true if it is multiplied throughout by e tB '. Thus, the whole principle 
is delivered by requiring that just the imaginary part is satisfied, but with C, p, o"o, P o 
replaced by Ce t9 , pe td aoe t6 , p^e for any 9. Henceforth, just the imaginary part of (13 .ip will 
be considered, while bearing in mind that C, p, <r and Po can be changed as just indicated. 
For the validity of the reasoning to follow, 9 should be chosen such that the imaginary parts 
of Ce t9 , and — pe t9 are positive definite. If p is real, as is typically the case, and C" is strictly 
positive definite, such a replacement allows the ensuing analysis to apply since (aside from 
in the variational principles (I3.24p and ( I3.26P ) we assume that both C" and — p" are strictly 
positive definite. 

The functional expressed in the imaginary part of (13. ip is saddle-shaped. An alternative 
functional is obtained by performing a partial duality transformation which involves a' and 
p" in place of e" and uu" . To clarify what is going on, note first that the time- reduced 
constitutive relations a = Ce, p = iupu imply that 
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Now following Cherkaev and Gibiansky (1994), elementary calculation gives 



He' e"\ 



C" a \ ( e' 
C -C" 



inf Ve" + |(e' o*)C[" 



a 



(3.3) 



where 

C" + C'(C")- 1 C -C'iC")- 1 

-{c")- x c (cy 1 

is positive definite, as can be seen from the positivity of the associated quadratic form 



C 



(3.4) 



e' a')C 



a 



e'C'e' + (CV - a'){C")-\C'e' - a'). 



(3.5) 



The infimum is attained when a' satisfies (I3.2p i . or equivalently a = Ce. Similarly, 



i ( urn' uiu" , , , „ 
; — p p 



p" —p' \ ( uu' 



^, 7/ -inf j-c^V + i(un/ v")V[^, \ \. (3.6) 



V 



(3.7) 



where 

V' + pVTV -p'ip'T 1 
-{ P ")- l p' (p'T 1 

is positive definite. Attainment of this infimum is consistent with ( 13. 2p ? or equivalently 
p = iupu. 

The imaginary part of (13.1 p can now be written 



inf inf sup 

u' a',p" U n j n 



e! e" 



-cr n 



a' -a' 



uu' uu" 



V" ~ Po 



(3.1 



Now 



sup f {e"(a'-a' )-uu"(p"-p'^)}dx = [ u^(a'-a' )ndS = [ {e'^(a'-a' )-uu'^p"-p'^} dx 

u" JQ JdQ JQ 

(3.9) 

so long as 

divcx' + f'= -up" in tt (3.10) 
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and a' satisfies any given traction conditions on dQ. Otherwise, the supremum is infinite. 
Therefore, the variational principle fl3.8j) reduces to 



inf inf 



e' e" 



-CTn 



a' -a' 



P Po 



+ i(e' a')C[ , )+h{um' p")V 



a 



UJU 

p" 



dx, 



(3.11) 



where p" is now given by (I3.10p and the infimum is subject to any prescribed traction 
boundary conditions on a' and any prescribed displacement boundary conditions on v! . 
Since the products e' W and u'qj/q do not affect the minimizer the variational principle further 
simplifies to 



inf inf 



a 



'Po uul. 



UJU 

p" 



+ i(e' a')C[ p")vH f 



dx. 



(3.12) 



Notice that the constitutive relations o = Ce and p = iuopu imply the new constitutive 
relations 



a 

-e" 



C 



a' 



P 

uu" 



V 



um 
P" 



(3.13) 



The fields appearing on the left sides of equations (13.131) are linked by the differential con- 
straints 

diva" + /" = tup', e" = \{Wu" + (Vu") T }, (3.14) 
while the fields on the right sides are linked by 



diva' + /' 



-cop", e' = §{W + (Vu'f}- 



(3.15) 



Following Milton, Seppecher, and Bouchitte (2009), we can rewrite (13 . 13f) as a single equa- 
tion G = CF where 



C 



C 
V 



G 



( °" \ 

-e" 

P' 

\ uu" J 



( e> \ 



a 



uu 
\P" J 



(3.16) 



and G satisfies the differential constraints (I3.14p while F satisfies the differential constraints 
(I3.15p . and C is positive definite. If we regard (I3.13p -f l3.15p as a primal problem then the 
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variational principle for it is 



inf I (—GqF + \F T CF) dx, 
F Jn 



(3.17) 



which is equivalent to (I3.12p . To check this variational principle directly, and thus provide 
an alternative derivation of fl 3 . X 2 1) . we compute the first order variation in the integral in 
(EH]): 



(G-G ) T 5Fdx 



[5u'(a" 



On )n 



(u^ - u")(5a')n] dS. 



(3.18) 



'n Jan 
This vanishes under the stated boundary conditions on the fields. It also clearly vanishes 
(and produces a valid variational principle) under more general boundary conditions. For 
example, although difficult to realize physically, one could fix the real part of the displacement 
and the real part of the traction at the boundary, so that 5u' = and (5<j')n = 0. Then 
choosing u' ' = 0, <r' ' = and up' Q = f" the variational principle ( I3.12p reduces to that of 



Milton, Seppecher, and Bouchitte (2009). Alternatively, by having no constraints on 8u' and 



(Sa')n one obtains a variational principle in which the trial fields are not required to satisfy 



any boundary condition, and at the minimum a"n = a'^n and u" 



u'q on dtt. 



If one prefers one can rewrite the variational principle in a form where only the surface values 
of u'q and a'^n appear. Using the differential constraints on u'q, e' Q \ (Jq and p' Q and integrating 
by parts, A3. 12[) is equivalent to 



inf inf 



-fu + He' a')C(0+i(um' p")vH, 



dx 



+ / (—uVo'n + u' Wn) dS. 
Jan 



(3.19) 



When / = the minimum of the variational principle (I3.12p . or equivalently (13. IT)) is 

/ \{G - 2G ) T Fdx = | f [u'(o" - 2o£)n + (2< - u")a'n] dS, (3.20) 
Jn Jan 

and thus depends only on the surface displacements u and u'q and surface tractions an 
and (TqU. Thus for any given choice of trial fields u' and & meeting the chosen boundary 
conditions, which imply that 

(yf - u')(a" - a r o)n = 0, « - u")(a! - a')n = 0, onffi, (3.21) 



8 



equations (13. 13|) and ( I3.2(jp produce the inequality 



dx 



> \ [ {(2u - u')a"n + u"(a' - 2q[)n] dS, (3.22) 
Jdn 



in which (I3.2ip has been used to eliminate a ' and Uq. Since <Tq and u' ' do not appear in this 



equation we may as well choose o'ln 



a"n and u'q 



u" on dVt in which case we are free 



to choose any yf and &. The inequality (j3.22j) bounds the possible surface fields u and cm, 
i.e. it bounds the Dirichlet to Neumann map associated with Q. Alternatively if the surface 
fields are known from measurements, then, for each set of trial fields, the inequality provides 
bounds on the moduli C and p inside Q and could potentially be used for tomography. 

In the limit as p" — > the last term in the variational principle ( 13. 12ft will remain finite if 
and only if 

p"/(M, (3-23) 



u 



and in this limit the variational principle becomes 



inf 



-<fA 



a 



UJUr 



oou 

p" 



a 



a 



dx, 



(3.24) 



where now p" and u' are given by (13.101) and (I3.23j) . 

These variational principles are particularly useful if a multiphase medium occupies Q and 
one of the phases, occupying a subregion \1/ is lossless, with C" = and p" = 0. Then 
assuming C and p' are positive in \1/ one does not have the flexibility to replace C and 
p by Ce td and pe ld since then the imaginary parts Ce ld and pe %e would lose their positive 
semi-definiteness within ^ . In the limit as C" — > in we see from ( 13. 5 j) that the last term 
in ( I3.24p will remain finite if and only if 



a 



C'e' in 



(3.25) 



In this limit the variational principle becomes 



inf 



-(Jo. 



.PO WM 



UJU 

p" 



dx 



a')C 



dx 



(3.26) 

in which p" and v! are given by ( 13. 1Q|) and (I3.23|) , and the infimum is over fields a' such that 
v! is an exact solution of the elastodynamic equations in i.e. ( 13 . 15j) . (I3.23p . and ( I3.25p . 
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4 Related elastodynamic principle of Hashin— Shtrikman 
type 

Recall that the constitutive relations (I2.2p are equivalent, in the case of time-harmonic fields, 
to those given in (13.131) . Now introduce a "comparison" medium, characterized by matrices 
Co, Vq (not necessarily obtained from complex tensors corresponding to C, p), introduce 
"polarizations" t", rj", tt' and z/', and define 

■>) <•)•(*•)• <"> 

These relations are consistent with the desired constitutive relations (I3.13P if the polarizations 
are chosen so that 

t" \ , , , / ( ' \ ( n' \ , m _ v / urn' 



(C-Co) , " =(V-V ) - ■ (4-2) 



-rj j \ o ) \ —v ) \ p 

With this preparation, note that 



+ (e> a')( yj-Hr" ~v" ) (C - C^ 1 ( ) . (4.3) 



Similarly, 



oju' \ , , , ^ I urn' 



+ (-m' ^l-KV -^')(P-^ ) )- 1 ( ^ ).(U) 



Equality is achieved when relations (14. 2 p are satisfied. 

The Hashin-Shtrikman variant of the minimum principle (13. lip is now 



inf inf 



o 



(r"-a ' e '-7/') l> 



a 

+ cr')Co( O-Ht" -?/' ) (C - Co)" 1 ( 
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urn 
p" 



i(tt' -i/)(7>-7> 



7T 



dx. (4.5) 



This produces an approximate solution for any given choice of the "polarizations" t", it', t]" 
and z/, and the exact solution if these are chosen to make the functional stationary. When 
C and V are chosen so C > C and Vo > V then the ~ signs in (14. 3 p and (I4.4p can both 
be replaced by < signs and then one has a minimum principle which involves taking the 
infimum over the "polarizations" r", tt', r)" and u', as well as an infimum over the fields 
u' and a'. If on the other hand Co and Vo are chosen so Co < C and Vo < V then the ~ 
signs can both be replaced by > signs and one obtains a saddle principle which involves 
taking the supremum over the "polarizations" and an infimum over the fields u' and a'. In 
either case the infimum over the fields v! and a' can be computed in terms of the relevant 
Green's function and one is left with a minimum principle or a maximum principle. When 
the body is infinite in extent the Green's function can be explicitly computed as shown in 
the appendix. 

In the notation of (I3.16P and (I3.17p . (14. 5 p becomes 



inf 

F 



[ [(T - Go) T F + \F T CoF - \T T {C - Cq^T] dx, 
Jn 



where 



T 



( t» \ 

-rj" 

7T' 



(4.6) 



(4.7) 



and in this form the principle is easily generalized to acoustics and electromagnetism, as will 
become clear in sections 5 and 6. It is perhaps worth remarking that there is no need for C Q 
to be restricted to have the block- diagonal structure displayed by C in (13 . 1 6j) i : this, however, 
is not developed further at present. 

Developing (14 .6p a little further, let F solve the variational problem for the comparison 
medium, i.e. in the case that T = 0, and write 



F = F + F'. 

Any admissible F' must satisfy homogeneous boundary conditions and hence 

F' T (-G + C F ) dx = 0. 



(4- 



(4.9) 
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It follows that (14. 6 p can be written 

inf / (-G^F + ±F T £ F + T T (F + F') + \F lT C F' - |T T (£ - C^T) dx. (4.10) 

Formally, if the field F' that attains the infimum is 

F' = -H T, (4.11) 

then f l4.10p becomes 

/ (-G^F + |F T £ F + T T F - \T T [(£ - Co)- 1 + Hq] 7 1 ) dx (4.12) 

and this is stationary with respect to T when 

[(£-£ )- 1 +Uo]T = F . (4.13) 
The operator "H relevant to an infinite medium is given explicitly in the appendix. 

5 Variational principles for acoustics 

For acoustics at fixed frequency in the presence of a body force / the fields satisfy the 
differential constraints 

-VP + / = zwp, h = V-v, (5.1) 

where P is the pressure, p the momentum, v the velocity, and h is defined to be the divergence 
of the velocity. They also satisfy the constitutive relations 

h = —iukP, v = rp, (5-2) 

where k = 1/k is the compressibility, the inverse of the bulk modulus k, which is complex 
with k" < 0, and r = is the inverse density (which we allow to be matrix valued and 
complex with r" > 0). Introducing the positive definite matrices 

/ k" + (k'f/k" -k'/k" \ fr" + r '(r")~V - r '(r")- 1 \ 

^ ~ I W IU" 1 IUII ' — I f*.l'\-l^l f„H\-l ' (5.o) 



-k'/k" 1/k" J 1 \ -{r")- l r' (r 

the constitutive relations take the new form 

h ' (-j')-nK), (5.4) 



-uP" ) \ h" J ' V P' J \ v 
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which can be expressed as the single equation G = CF with 



C 



K 
K 



G 



( hi \ 
-ujP" 
—ujv' 

\ up' J 



(-ujP'\ 
h" 
up" 
\ ujv" ) 



(5.5) 



where the fields appearing in G are linked by the differential constraints 

U = V • v', -VP" + /" = cop', 
while the fields appearing in F are linked by 

- VP' + /' = -up", h" = V • v". 



(5.6) 



(5.7) 



(The second equation in (15 .4p was written with the real components of the fields on the 
left, rather than on the right, to ensure that the differential constraints couple fields on the 
same sides of (I5.4p .) Based on (I3.17p . we assert that the form of the variational principle for 
acoustics is 



inf inf 

v" P' 



n L 



+ U-uP' h")K 



+ ( v' Q -p' ) 



p" 



-ujP' 
h" 



+ \u 2 {p" v")K 



P" 



dx, (5.8) 



in which h" and p" are given by (15.71) . and h' , Pq, v' , and p' satisfy the differential constraints 
(15.61) . To check this we see that the first order variation 



h'-h' cjPH 



uP" 



-u5P' 
Sh" 



+ u 2 (v' -v' p' -p' ) 



dp" 

Sv" 



dx 



u / [5P'(v' - v')n + (P^ - P"){5v")n] dS, 
Jan 



(5.9) 



only depends on boundary terms. (Note the need for the factor oo 2 in (I5.9P which explains 
the need for the additional factors of uo in the last two components of G and F in (15.51) ). 

The condition that the boundary integral vanishes for all 8P' and Sv" that are permitted 
according to which components of P' and v"n are prescribed on dfl, forces the complementary 
components of (v' — v')n and Pg' — P" to be zero, thus fixing these components of v'n and 
P" . In the particular case when / = and all components of P' and v"n are prescribed at 
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the boundary, then by choosing Pq = 0, p' = 0, and v' = 0, these variational principles 
reduce to those of Milton, Seppecher, and Bouchitte (2009). 

If one prefers one can rewrite (I5.8P in a form where only the surface values of Pq and v' n 
appear. Using integration by parts the variational principle is equivalent to 



inf inf 

v" P> 



-uP' 
h" 



+ \u 2 {p" v")K 



dx 



—uf v + i ( —uP' h")K 

+ u I {P'v' n + Po'v"n) dS. (5.10) 
Jan 

When / = the minimum value of (I5.8P or (I5.10p only depends on the boundary fields, and 
is 

/ [P'(2v' -v')n + (2P^-P")v"n]dS. (5.11) 
Jan 



Thus knowledge of this boundary data will, for any given choice of trial fields v" and P', 
give constraints 



-ujP' 



uP! h»)K[ It )+^ 2 (p" v")K(Z 



dx 



> 



\u f [{P'-2P!)v'n + P"(v"-2v")n) dS (5.12) 
Jan 



on the possible values of the moduli k and p inside Q, and thus may be useful for tomography. 
In the limit as r" — > 0, i.e. as p" — > the last term in (15. 8p will be minimized when 



v" = r'p", 



and is infinite otherwise. In this limit the variational principle reduces to 



inf 



-luP' 



-h' coP^i - j+^(v' -p^i^j+hi-uP' h")K 
where h" and p" are given by (15.71) . while v" is given by (I5.13p . 



-uP' 
h" 



(5.13) 



dx, 
(5.14) 



6 Variational principles for electrodynamics 

For electromagnetism it is more conventional to have fields which are multiplied by e~ tult , 
rather than e luJt , so that the permittivity tensor e and permeability tensor \i have positive 
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definite imaginary parts. We now adopt this convention. Then the fields in Maxwell's 
equations at fixed frequency u satisfy the differential constraints 



Vxfi = icoB, V x H = -icuD + j, 
and the constitutive equations 

D = eE, H = mB, 
where m = jjT 1 . Introducing the positive definite matrices 



£ 



e'(e"Y x e' 



-e'ie'T 1 



M 



{s")- 1 

the constitutive relations take the new form 

E' 



m" + m'(m") l m' 



-m!(m") 1 



(m 



rr\-l 



D" 
—E" 



£ 



H' 
-B' 



M 



D' ' 

which can be expressed as the single equation G = CF with 



B" 
H" 



C 



£ 
M 



G 



( D" \ 




(E'\ 


—E" 




D' 


H' 


, F = 


B" 


\-B'J 




\H") 



(6.1) 
(6.2) 

(6.3) 
(6.4) 



(6.5) 



where the fields appearing in G are linked by the differential constraints 

VxH'-j' = luD", V x E" = luB', 
while the fields appearing in F are linked by 

V x E' = -ojB", V x H" - j" = —luD'. 
Based on (I3.17P we state that the form of the variational principle for electromagnetism is 

-(-H' B' ) 



(6.6) 
(6.7) 



inf inf 

E' H" 



E' 
D' 



B" 
H" 



+ D')S[ %)+HB" H»)m(%'„ 



dx. 



(6.1 



where D' and B" are given by (j6.7j) . while Dq, Eq, H' q and B' satisfy the differential con- 
straints ( 16. 6L To check this we use the identities 

V • (E' x H') = H' ■ (V x E') - E' ■ (V x H') = -uo{H'B" + E'D") - j'E', 

V • (E" x H") = H" ■ (V x E") - E" ■ (V x H") = u(H"B' + E"D') - j"E", (6.9) 
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which hold for any fields satisfying the differential constraints (16. 6p and (16. 7p . Then we see 
that the first order variation can be expressed in terms of boundary fields: 



D" — Dq Eg-E" 



5E' 



5D 1 



, + (H'-H> B>-B'[ 



SB" 
5H" 



dx, 



+V ■((#'- H' ) x 8E') +V-((E'o- E") x 5H") 

- [ [((#' - H' ) x 8E')n + (« - £") x <fif") n ] dS. 
w Jen 



dx. 



(6.10) 



Note that this surface integral only depends on the tangential components of the fields 
H' — Hq, SE', Eq — E", and 5H" . The condition that the boundary integral vanishes for all 
8E' and 5H" that are permitted according to which tangential components of E' and H" 
are prescribed, forces the complementary tangential components of H' — H' and Eq — E" 
to be zero, thus fixing these tangential components of H' and E" . Again these variational 
principles reduce to those of Milton, Seppecher, and Bouchitte (2009) when the tangential 
components of E' and H" are fully prescribed. 

If one prefers one can rewrite (16. 8 ft in a form where only the surface values of H and E' ' 
appear. Using integration by parts the variational principle is equivalent to 



inf inf 

E' H" 



E' 

D' 



UB" H")M 



B" 
H" 



dx 



j r E'/u + i(E' D')8 

+ / [-(H' x E')n + « x H")n] dS. 
Jan 

When j = the minimum value of ( 16. 8ft or (16.111) is 

/ [((H' - 2H ) x E')n + ((2E' ' - E") x H")n] dS. 
Consequently for any given choice of trial fields E' and we have the inequality 



(6.11) 



(6.12) 



> 



dx 



/ l(H' x (2g - E'))n + (E" x (H" - 2H!'))n] dS, (6.13) 
2^ Jan 

which could be useful for tomography if we have measurements of the boundary values of 
the fields. 
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In the limit as //' — > 0, i.e. as m" — > 0, the last term in the 

H" = m'B" 



will be minimized when 

(6.14) 



and is infinite otherwise. In this limit the variational principle reduces to 

E' 



inf 



n L 



E' 
D' 



B" 
H" 



\{E' D>)£ 



D' 



where D' and B" are given by (16. 7p . while H" is given by (I6.14p . 



dx, 



(6.15) 



Appendix: Green's function for infinite comparison medium 

The Hashin-Shtrikman variational functional is rendered stationary (with respect to the 
fields e', a', urn', p") when e" is the strain field associated with the displacement u" and the 
equation of motion 



div 



:■)•(> 

is satisfied, together with the relations (14. ip . 
Motivated by the form of relations (13. 2p . define 



CO 



D 2 D l 



-D* 



and Q 



p 

-p" 



Q2 Qi 
QI -Qs 



(A.1) 



(A.2) 



in which D2, -D3, Q2 and Q3 are symmetric real matrices and, in terms of these and the real 
matrices D\ and Qi, 



C ° ~ 1 -D^DJ ) and Vo ~ ~ { -QtQl Q3- 1 



The relations (14. ip may then be given in the equivalent form 



a" - t" 



a 



1 ' \ / p' — n' 



1 e" - 77" ' ' 



-p" 



Q. 



1 uu" + v' 



(A.3) 



(A.4) 



The quadratic forms associated with Co and Vq are positive-definite if D2, -D3, —Q2 and —Q3 
are positive-definite, as will henceforth be assumed. 

The equation of motion (lA.ip now gives 

div(£> W) +u 2 Q U + f = 0, (A.5) 
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where 



U=[ i) and / = ( f" + W «< "*''). (A.6) 



The required Green's function Q is that associated with equation flA.5[) . It can be con- 
structed by a method introduced for the "ordinary" elastodynamic Green's function by 



Willis (1980). First, select f(x) = foS(x), with f constant, and note the plane-wave decom- 



position of the three-dimensional delta-function 



6 ( x ) = ~SZa / " X ) dS - ( A - 7 ) 

iei=i 



8tt 2 



This motivates first solving (A.6), with / replaced by /o5(£ -x), where |£|=1. The solution, 
say, depends on x only in the combination s = £ • x; thus, it satisfies the system of ODE's 



<Pu* 

ds 

where 



4>^T + " 2 QoU^ + f S{s) = 0, (A.8) 



4,= (£n£ C^)- (A.9) 



.£' />;£ 

The system can be diagonalized by reference to the eigenvalue problem 

C U N = c 2 N Q U N (N = 1,2,- • -6). (A.10) 

Although the matrices are real and symmetric, they are not positive-definite and the eigen- 
values may be real or they may occur in complex conjugate pairs. Suppose first that (IA.10I) 
is satisfied for some real c N and (hence real) Un, having the form shown in equation (13)i: 
say Un = (u N ,u" N ) T , where u' N and u" N are 3- vectors. The 3- vector u" N can be eliminated 
from (1Q0]| . to give 

{[(£ A£) - c%q 2 ] + [(cao - <? n Qi][{zd 3 £) - 4g 3 ]^[(^i0 - <? n Qi] t }u'n = o (A.n) 

and hence 



[v. 



NJ 



{[(^0-4Q2]+[(eA0-4Qi][(^30-4g 3 ]- 1 [(^i0-4Qi] T }"'iv = o. (A.12) 



Recalling the definiteness of D 2 , D 3 , —Q 2 and — Q%, it follows that the expression on the left 
side of (1A.12|) cannot be zero if c N is positive. Hence, any real eigenvalue must be negative. 
Now write 

cl = (c' N + id' N f , N = !,■■■&. (A.13) 
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Since c" N cannot be zero, it is possible to ensure that c" N > for all N. If, furthermore, c 2 N 
is complex, then c 2 M = c 2 N and Um - 
eigenvectors are normalized so that 



is complex, then c 2 M = c 2 N and Um = C/jv f° r some M, and correspondingly cm = —cn- The 



Setting 



it follows that 



UJjQqUm — 8nm- 



M=l 



(i | • 11 \2 
(C N -+- IC N ) 



d 2 <p N 
ds 2 



+ co 2 4> N + U^f S{s) 



0. 



(A. 14) 



(A.15) 



(A.16) 



The solution that decays to zero as \s\ — > oo is 

(p N (s) 



U£f exp 



2iuj(c' N + ic" N ) 



Having now determined U^(s), the solution U(x) of ( 1A.5I) . with f(x) 
making the superposition implied by ( 1A.7I) : 



where 
Go(x) 



E 



,AT=1 



u N u% 



2| ^ [ - ( c' n + ic 



U{x) = Go(x)f , 



6{t ■ x) 



(A.17) 
fo6(x), follows by 

(A.18) 



^) 2 



exp 



AT; 



— IU)\£ ■ x\ 

<.Cn ~f~ ic 



N/ 



dS 



(A.19) 

If, for any N, the eigenvalue and eigenvector are real, then c' N = and the contribution to 
the sum is real. If the eigenvalue is complex, it follows from the properties established above 
that some other M contributes the complex conjugate expression. Thus, Qq is real. 

Now for some general f(x), and in particular the f(x) given by ( 1A.6I) . the solution is 



U(x) 



Go(x-y)f(y) dy 



which gives us u' and u" and hence e' and e", and from (IA.4I) we get 

a' = Die' + D 3 ( V " -e"), p" = Q%uu' - Q 3 (uju" + v') 



(A.20) 



(A.2i; 
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Ho 



The functions u', a', e', and p" thus obtained in terms of the polarization fields are then 
substituted back in (14. 5p . The resulting functional is a minimization (or maximization) 
principle for the polarization fields when Co and Vq are chosen so both Co — C and Vo — V 
are positive definite (respectively, negative definite). 

This appendix is concluded by recording the explicit form of the operator "H introduced in 
(14. lip . For this purpose, it is helpful to define 

S i]k = \{5 ik d j + 5 jk d l ) (A.22) 

so that the strain e associated with displacement u is e = Su, and also to define 

^kji = £ijki (A. 23) 

but with the partial derivatives acting backwards, with respect to the immediately preceding 
variable, so that, for example, 

(G£t) _i{^ + ^}. (A , 4) 

We also write 

— G3 ) ' (A. 25) 

In terms of the representation (14. lip employing the polarization T as defined in (14 .7p . it 
follows from fTA~5]) . fTA~20]) and flA~2T]) that 

£G 2 tf S(G 2 £*>D 1 -G 1 SW 3 ) u>£G 2 w£(G 2 Qi - G1Q3) 

{Df £G 2 - D 3 £GJ )£t |d 3 + Df£{G 2 £ t£>j - Gi£t£> 3 ) w(Df £C? 2 - D 3 £Gf ) £(G 2 Qi - GiQ 3 ) 

-D 3 £(Gj£tDi + G 3 £tL> 3 )} -D 3 £(GfQi + G3Q3)} 

uG 2 £t uiiG^D! - Gi£+Z) 3 ) w 2 G 2 ^ 2 (G 2 Qi - G1Q3) 

w(QfG 2 -Q 3 Gf)£t w|Qf(G 2 £tD 1 -Gi£tD 3 ) w 2 (Qf G 2 - Q 3 Gf ) -{q 3 - u> 2 Qf (G 2 Qi - GiQ 3 ) 
-Q 3 (Gf f+Di + G 3 f tD 3 )} WQ 3 (Gf Qi + G 3 Q 3 )} 

(A.26) 

The symmetry of Ho is consistent with its association with the variational formulation for 
the comparison medium. 
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